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We have performed experiments using a 3D-Bose-Einstein condensate of sodium atoms in a ID 
optical lattice to explore some unusual properties of band-structure. In particular, we investigate 
the loading of a condensate into a moving lattice and find non-intuitive behavior. We also revisit 
the behavior of atoms, prepared in a single quasimomentum state, in an accelerating lattice. We 
generalize this study to a cloud whose atoms have a large quasimomentum spread, and show that 
the cloud behaves differently from atoms in a single Bloch state. Finally, we compare our findings 
with recent experiments performed with fermions in an optical lattice. 
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An optical lattice is a practically perfect periodic po- 
tential for atoms, produced by the interference of two or 
more laser beams. An atomic-gas Bose-Einstein conden- 
sate (BEC)ll is a coherent source of matter waves, 
a collection of atoms, all in the same state, with an ex- 
tremely narrow momentum spread. Putting such atoms 
into such a potential provides an opportunity for explor- 
ing a quantum system with many similarities to electrons 
in a solid state crystal but with unprecedented control 
over both the lattice and the particles. In particular we 
can easily control the velocity and acceleration of the lat- 
tice as well as its strength, making it a variable "quan- 
tum conveyor belt" . This allows us to explore situations 
that are difficult or impossible to achieve in solid state 
systems. The results are often remarkable and counter- 
intuitive. For example atoms that are being carried along 
by a moving optical lattice are left stationary when the 
still-moving lattice is turned off, in apparent violation of 
the law of inertia. 

A few experiments have studied quantum degenerate 
atoms in moving optical lattices ^,^,^,^]. Bragg diffrac- 
tion of a Bose condensate is a special case of quantum 
degenerate atoms in a moving lattice 0. Here, using 
a Bose-Einstein condensate and a moving lattice, we 
achieve full control over the system, in particular its ini- 
tial quasimomentum and band index as well as its subse- 
quent evolution. We also show the difference in behavior 
when the atom sample has a large spread of quasimo- 
menta, as compared with the narrow quasimomentum 
distribution of a coherent BEC. 

Our lattice is one-dimensional along the x axis, pro- 
duced by the interference of two counter-propagating 
laser beams, each of wave- vector k = 27r/A (A ^ 589 
nm is the wavelength of the laser beams). This results 
in a sinusoidal potential, V s\v? kx, with a spatial period 
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A/2. 

We will use Bloch theory, emphasizing the single par- 
ticle character of the problem. An overview of Bloch 
theory, as it applies to this one dimensional system, is 
supplied in reference Q . Briefly, the wave function of the 
atoms in the lattice can be decomposed into the Bloch 
eigenstates Un,q{x)e'^^^ characterized by a band index n 
and a quasimomentum defined in the rest frame of the 
lattice. The eigenenergies of the system, En{q)^ as well as 
the eigenstates are periodic in q with a periodicity 2hk, 
the reciprocal lattice vector of the lattice. A wave packet 
in band n with quasimomentum distribution centered at 
has a group velocity Vg = dEn{q)/dq. Fig. Q shows 
the band structure in the repeated-zone scheme 8], for 
a lattice with a depth V = AE^ {E^ is the single-photon 
recoil energy given by E^ = h^k'^/2M and is related to 
the recoil velocity by = M being the mass 

of an atom) . Note that for convenience the band-energies 
En are offset such that they coincide at large band in- 
dex with the free parabolae; this shows more clearly the 
avoided crossings between free particle states due to the 
laser-induced coupling. These avoided crossings create 
the band gaps that separate energy bands with different 
indices n. 



I. EXPERIMENTAL SETUP 

The experimental setup has been described previ- 
ously 0. An almost pure Bose-Einstein condensate (no 
discernable thermal component) of about 2 x 10^ sodium 
atoms is prepared in a triaxial Time Orbiting Potential 
(TOP) trap We adiabatically expand the conden- 

sate by lowering the mean trapping frequency 20] from 
200 Hz to a value ranging from 100 Hz to 19 Hz. This 
reduces the atom-atom interaction, the strength of which 
is given by the chemical potential /i = ^^^^^ ^^], n be- 
ing the density at the center of the cloud and a ~ 2.8 
nm the scattering length. During the expansion, the cal- 
culated Thomas-Fermi diameter, 2i?, of the condensate 
along the lattice direction increases from 18 /im up to 
values ranging from 24 /im to 48 /im. The wave-function 
of each atom thus covers more than a 100 lattice sites and 
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FIG. 1: Band structure for a 4^r deep lattice in the repeated- 
zone scheme. The dotted hues represent the free particle 
parabolae to which the bands adiabatically connect as V ^ 0. 
The region in dark grey corresponds to the first Brillouin zone. 
The region in light grey corresponds to the second Brillouin 
zone. 



is an excellent approximation of a Bloch state. The rms 
width of the momentum distribution of the atoms in the 
condensate along the axis of the lattice is 
Therefore the rms width of the quasimomentum distri- 
bution of each atom is Aq ~ (X/AR) fik < O.Olfi/c. 

To form the lattice we use two counter-propagating 
laser beams perpendicular to gravity. Each has a power 
of up to 10 mW and is detuned either 200-350 GHz to 
the blue of the sodium D2 transition (experiments of sec- 
tions HH and HVj) or 700 GHz to the red of the Dl tran- 
sition (last section). They are focused to a beam waist 
of about 200 jiia FWHM, leading to a calculated sponta- 
neous emission rate < 30 s~^, negligible during the time 
of the experiments. The lattice depth, measured by ob- 
serving the Bragg diffraction j^, is up to 13£^r- We use 
acousto-optic modulators to independently control the 
frequencies and intensities of the beams. The unmodu- 
lated intensity is kept constant to within 5% by active 
stabilization. A frequency difference 5 between the two 
beams produces a "moving standing wave" of velocity 
V = 5/ 2k. Numerically, a difference of J/27r = 100 kHz 
corresponds to a lattice velocity of one recoil velocity, 
Vr = hk/M ^ 3 cm/s. 

The cloud's momentum is analyzed using time-of- 
flight. The time-of- flight period, typically a few millisec- 
onds, converts the initial momentum distribution into 
a position distribution, which we determine using near- 
resonance absorption imaging along an axis perpendicu- 
lar to the axis of the lattice. 



II. DRAGGING A CONDENSATE IN A 
MOVING LATTICE 

In a first set of experiments, we begin with a BEG 
in a magnetic trap with a 19 Hz mean frequency. This 
weak trap makes the interactions between atoms almost 
negligible on the time scale of the experiment, i.e. h/uAs 
generally longer than the duration of the experiment [2l| . 
After turning off the magnetic trap, we adiabatically ap- 
ply a moving lattice with a final depth of AE^. The 
turn-on time of the lattice intensity is 200 /is, an interval 
chosen to ensure adiabaticity with respect to band ex- 
citation j22j (see section nil)) . The fixed velocity of the 
lattice, V, is between and about 3vr. In the lattice frame 
the atoms have a quasimomentum q = —Mv. Because 
the width of the quasimomentum distribution is very nar- 
row, this procedure produces a good approximation of a 
single Bloch state with a freely chosen q. 

Atoms loaded in this way are dragged along with the 
moving lattice. In the limit that the lattice is very deep 
so that the bands are flat (i.e. = 0), the group 

velocity with respect to the lattice, Vg, is and the atoms 
are dragged in the lab frame at the velocity of the lattice. 
For finite depth lattices the dragging velocity in the lab 
frame is v -\-Vg. (Note that for > 0, '^g < so that this 
dragging velocity in the lab frame v -\- Vg < v.) 

In order to experimentally measure the dragging ve- 
locity we suddenly (on the order of 200 ns) turn off the 
moving lattice, projecting the Bloch state onto the basis 
of free-particle momentum eigenstates while preserving 
the momentum distribution. Figure (SJi shows the lat- 
tice depth as a function of time. Images of the resulting 
diffraction pattern for various lattice velocities are pre- 
sented in figure [2l3. The average velocity seen from the 
diffraction pattern (the weighted average of the veloci- 
ties of the individual diffraction components) increases 
with the lattice velocity through the first Brillouin zone. 
In fact for this rather flat band the dragging velocity is 
roughly equal to the lattice velocity. (The details for 
higher velocities are discussed in the following section.) 

An alternate method to study the atomic momentum is 
to release the condensate adiabatically (~ 200 /is) rather 
than suddenly, thus avoiding diffraction. Figure|3Ji shows 
the lattice intensity time sequence for this method. The 
corresponding images for various lattice velocities appear 
in figure int. These pictures show that (apart from when 
the lattice velocity is very close to an integer multiple 
of Vr, a situation discussed in section ITTH) the atoms are 
back at rest in the laboratory frame, despite the fact that 
the lattice is still moving during the ramping down of its 
intensity. This is true even in the first Brillouin zone 
where the lattice drags the condensate at roughly the 
lattice velocity. This result is especially surprising when 
one considers that atoms moving with the lattice return 
to zero velocity as if they had no inertia. One might also 
ask how do the dragged atoms "know" that they should 
be at rest when the lattice is turned off. One way of 
understanding this is to note that the lattice turns on 
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FIG. 2: (Color) Dragging of atoms in a moving lattice fol- 
lowed by a sudden turn-off of the lattice. Fig. presents the 
time sequence. Fig. 13) shows the absorption image of the 
cloud after a 1.5 ms time of flight following a sudden turn off 
of the lattice for different lattice velocities v, related to the 
quasimomentum by ^ = —Mv. The numbers on the vertical 
axis refer to the atomic velocity in units of Vr. The average 
velocity of the atoms in the lab frame, deduced from fig.l^t, is 
shown in fig. 13: versus the velocity of the lattice. The initial 
velocity of the condensate in the magnetic trap fluctuates with 
an RMS value of O.OSi^r- The mean velocity, after suddenly 
turning off the lattice, thus exhibits the same fluctuations. 
The solid curve is the mean velocity of the atoms calculated 
from the band structure for a 4:Er deep lattice. 
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FIG. 3: (Color) Dragging of atoms in a moving lattice followed 
by an adiabatic turn-off of the lattice. Fig. EJt presents the 
time sequence. Fig. |3) shows the absorption image of the 
cloud after a 1.5 ms time of flight following the adiabatic turn- 
off of the lattice for different lattice velocities. The numbers 
on the vertical axis refer to the atomic velocity in units of Vr- 



adiabatically and turns off adiabatically along the same 
path. This must necessarily return an eigenstate of the 
hamiltonian to the same eigenstate. A more detailed 
explanation involving band-structure will be presented 
in the next section. 



III. ANALYSIS OF THE EXPERIMENTS 

All experiments described in this paper start with an 
adiabatic turn-on of a lattice moving at a velocity v. In 
the lattice frame, in the limit of a vanishingly small lattice 
depth, the atomic wavefunction of a momentum eigen- 
state has a phase gradient —Mv/fi corresponding to the 
velocity —v of the atoms with respect to the lattice. This 
free particle state is also a Bloch state with a quasimo- 
mentum q corresponding to a phase gradient q/?i^ so that 
q = —Mv. All changes in the lattice intensity preserve 
this quasimomentum (as can be seen by calculating that 
the matrix elements of the periodic potential between 
Bloch states of different q, are zero). When the lattice 
is fully turned on, the quasimomentum is still —Mv and 
if the turn-on has been adiabatic (so that no change in 
band index occurs), we end up in a single Bloch state. 

Referring to figure H we see that, when atoms are 
loaded adiabatically into the lattice with the quasimo- 
mentum in the first Brillouin zone, the free particle mo- 
mentum connects to the corresponding quasimomentum 
in the lowest, n = 0, band. For quasimomenta outside 
the first Brillouin zone, the free particle momenta con- 
nect to the corresponding quasimomenta in the appro- 
priate band. For example if the velocity of the lattice is 
1.5^'r, i.e. in the second Brillouin zone, the atoms will end 
up in the second, n = 1, band with a quasimomentum 
q = —1.5hk. There is thus a strict relation between the 
range of quasimomenta and the band index into which 
the atoms are loaded: if the quasimomentum is in the 
^th BriHouin zone, the atoms are loaded into band n — 1. 
On the other hand, if for example we wish to prepare 
the atoms in q = —1.5fik and n = 0, we would have to 
accelerate the lattice, as described in section Hvl 

The condition for adiabaticity with respect to band ex- 
citation during the loading has been detailed in ref. |3] : in 
order to avoid transitions from a given band to an adja- 
cent band, the rate of change of the lattice depth V must 
fulfill ^\{n,q\sm^kx\n ± l,q)\ <C AE^/h. AE is the 
energy difference between the given band and its nearest 
neighbor. When AE approaches (as is usually the case 
near a Brillouin zone boundary when V ^ 0) the process 
cannot be adiabatic. For (7 = 0,n = 0, AE > AE^^ and 
the natural time scale for adiabaticity with respect to 
band excitation is on the order of h/AEr. We emphasize 
that in the limit of V ^ there is a natural energy gap 
due to the periodicity of the lattice, AE ^ (except at 
the edge of the Brillouin zones). The existence of this 
non-zero energy gap when the lattice depth goes to zero 
is in contrast to, for example, a harmonic oscillator for 
which the spacing between energy levels does go to zero 
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as the strength of the potential vanishes. 

We now analyze in more detail the two methods for 
studying the momentum distribution described in the 
previous section. 

In the first method we turn off the lattice potential sud- 
denly, i.e. diabatically. This sudden turn-off leaves the 
atomic momentum distribution unchanged from what it 
was in the lattice. If the atoms are in a Bloch state, cor- 
responding to a single value of the wave-function as 
viewed in the rest frame of the lattice is a superposition 
of plane waves with momenta q + 2m?ik {m is an inte- 
ger) . The population- weighted average of the momentum 
components gives the mean momentum of the atoms in 
the lattice, which is Mv^ 8]. In the laboratory frame 
these momentum components are shifted by the velocity 
of the lattice and are observed as a diffraction pattern. 
The time-of-flight spatial distribution of these momen- 
tum components is analogous to the diffraction pattern 
of any wave from a periodic structure. The spacing be- 
tween the momentum components gives the reciprocal 
lattice vector, 2?ik in our case. This diffraction is char- 
acteristic of sudden turn-off (or on) of the lattice. 

Figure shows the measured dragging velocity in the 
lab frame as a function of the lattice velocity v = —q/M. 
Also shown is the calculated dragging velocity + v 

for a AE^ deep lattice. When the atoms are in the first 
Brillouin zone and in the n = band they are dragged 
along at close to the lattice velocity, because the n = 
band is nearly flat (see figure Q). The next, n = 1, band 
is much less flat and the atoms are not dragged at the 
lattice velocity except at the edge of the Brillouin zone 
where = vanishes. In the third Brillouin zone, 

the n = 2 band is so close to a free particle that there is 
almost no dragging and experimentally we do not even 
see good dragging near the zone boundary at Sv^ because 
the feature is too narrow. This behavior is rather intu- 
itive in that the lattice drags atoms effectively up to a 
velocity for which the atomic kinetic energy in the lattice 
frame is about equal to the lattice depth. Reference 4] 
reported similar results, measuring the dragging velocity 
using the displacement of the cloud rather than diffrac- 
tion. (Note that they plot the group velocity.) This 
dragging process is also discussed in |^. 

Now let us consider the rather counter-intuitive results 
obtained by adiabatically ramping off the lattice inten- 
sity. As noted earlier, turning off the lattice either adia- 
batically or non-adiabatically does not change the quasi- 
momentum distribution, although it may change the mo- 
mentum distribution. (This assumes that no other forces 
besides the lattice act on the atoms in the rest frame of 
the lattice. This assumption would be violated, for ex- 
ample, in the presence of interaction between the atoms 
or if the lattice were accelerated.) 

Consider a single Bloch state in the lattice, as is the 
case in the previous section. In contrast to the sudden 
turn-off method described above, the multiple momen- 
tum states q + 2m?ik coalesce into a single momentum 
component, whatever the depth of the lattice was. Look- 



ing at figure d we can see that any single Bloch state 
\n^q) will adiabatically connect to a single free particle 
parabola, unless there is a degeneracy and adiabaticity 
fails. For the specific experiment described in section ITTl 
where a lattice moving at a constant velocity is turned on 
and off, this parabola is always the one labelled Ofi/c. In 
this case the Bloch state produced is such that the single 
momentum component \s p = q = —Mv in the frame of 
the lattice. Transforming into the lab frame we find the 
velocity of the atoms to be zero, as observed. 

As an alternate explanation we recall posing the ques- 
tion "how do the dragged atoms 'know' that they should 
be at rest when the lattice is turned off?". We now can 
see that this information is stored in the phase gradient, 
or the quasimomentum, which does not change as the 
lattice is ramped on and off. We again emphasize that, 
in the absence of interactions, this phase information is 
preserved no matter how deep the lattice was or how fast 
the lattice was turned on and off. 

Let us now return to the failure of adiabaticity near the 
edge of the Brillouin zones. Referring to figure^ consider 
free atoms, stationary in the lab frame, but at the edge of 
a Brillouin zone in the lattice rest frame, for example at 
q = hk OT q = 2hk. At q = hk atoms will, as the lattice 
is turned on, be loaded into both bands n = and n = 1; 
at g = 2hk atoms will be loaded into n = 1 and n = 2. 
Upon turning off the lattice, the two populated states will 
each connect to two free-particle parabolae. For example 
dit q = hk atoms will be in both the Ohk and 2hk parabo- 
lae (at q = 2hk^ they will be in both Ohk and Ahk). In 
the lattice frame (with the lattice off) atoms at q = hk in 
the Ohk parabola are moving with a group velocity -\-Vr; 
atoms dit q = hk in the 2hk parabola are moving with 
a group velocity —v^. Transforming back into the lab 
frame, these atoms are moving at and —2vr respec- 
tively. Similarly at q = 2hk in the lattice frame the atoms 
are moving at -\-2vr and — 2vr, corresponding to and Av^ 
in the lab frame. This is exactly what is experimentally 
seen in figure IHId. (And it is exactly the same as first- 
and second-order Bragg diffraction 7]). The fraction of 
population in each momentum component depends on 
the details of the loading and the unloading. For higher 
bands the adiabaticity condition becomes easier to sat- 
isfy near a band edge. Even though the band gap at the 
level anti-crossing at the Brillouin zone edge gets smaller 
for larger band index, the energy difference AE between 
adjacent bands at a fixed distance in quasimomentum 
from the Brillouin zone edge, is larger for higher bands 
(see fig.^). This larger AE leads to greater adiabaticity 
for a given rate of change of the lattice depth at given 
distance in quasimomentum from the zone edge. This 
partially explains why so little population in non-zero 
momentum states is seen near the band edges for high 
velocities in fig. In addition the coupling between 
adjacent bands gets smaller for higher bands (because it 
represents a higher order process), as reflected by the nar- 
rowing of the band gap, and this smaller coupling further 
reduces the population of non-zero momentum states. 
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This method to analyze the quasimomentum distribu- 
tion by adiabaticahy ramping down the amphtude of the 
lattice is independent of the way this distribution has 
been created, and thus allows the analysis of complex 
quasimomentum distributions. In order to understand 
this point, we recall that there is a unique correspon- 
dence (except at the Brillouin zone boundary) between 
any given Bloch state \n^q) in the lattice and a momen- 
tum state in the lab frame when the lattice is adiabati- 
cahy turned off. For example two Bloch states with the 
same q {{) < q < Tik) in a lattice moving with a veloc- 
ity V, but in two adjacent bands, let's say n = and 
n = 1, will connect to momentum states q + Mv and 
q + Mv — 2hk respectively. In the same way two Bloch 
states with the same band index and two different g's 
will end up in two different momentum states. Suppose 
now we prepare a given quasimomentum distribution in 
the lattice frame, consisting of many q^s in many bands, 
and suppose we adiabaticahy ramp down the intensity of 
the lattice. If during that ramping-down time the quasi- 
momentum distribution does not significantly evolve (e.g. 
under the influence of interaction, or under acceleration), 
the adiabaticity ensures that the population in a given 
state of quasimomentum q in band n is conserved during 
the process. The quasimomentum distribution is thus 
mapped onto a momentum distribution in the labora- 
tory frame 111. This method, which has also been used 
in reference 12], then allows us to fully reconstruct the 
quasimomentum and band distribution. 

We will give other examples of such mappings in the 
next two sections. 



IV. ACCELERATION OF A CONDENSATE IN 
A SINGLE BLOCH STATE 

In this section, we revisit the behavior of atoms under 
acceleration of the lattice, already studied in 0,0,0, 
using the adiabatic ramp-down analysis described in the 
last section. For this particular experiment, we again 
decrease the mean oscillation frequency of the magnetic 
trap to 19 Hz before turning the trap off. Starting with 
the condensate at rest in the lab frame, we linearly turn 
on the stationary lattice intensity over 40 fis in order to 
ensure adiabaticity. The final depth for this experiment 
is y = ISEj-. All the atoms are now approximately in 
the state \n = O,^' = 0). We then accelerate the lat- 
tice for 400 /is up to a given velocity 'Uf, with a constant 
acceleration a < 800 m/s^. The quasimomentum q of 
the atoms evolves during the acceleration according to a 
lattice version of "Newton's law" q = —Ma 0- In the 
lattice frame this is equivalent to adding a linear poten- 
tial -Max. Provided that \Ma{l, q\x\0, q)\ <^ Ei - Eq 
|2^ . there is no transition between the first two bands 
and the atoms stay in the lowest band. This implies that 
the acceleration should be smaller than 4 x 10^ m/s^, a 
condition well satisfied in our experiment. This accel- 
eration allows us to produce any q in the lowest band. 
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FIG. 4: (Color) Acceleration of atoms in the ground state 
(n = 0) band starting from ^ = 0. The lattice is adiabaticahy 
raised up to ISEr, accelerated and then adiabaticahy turned 
off at constant velocity. Fig. 2Jl shows images of the conden- 
sate after time of flight for increasing final lattice velocities. 
Fig. 2J) shows the position of the center of mass of the cloud 
in the lab frame. The circles are the positions measured on 
fig-Eti whereas the crosses represent the position of the cloud 
minus the displacement due to the dragging of the lattice. It 
therefore gives the momentum of the atoms. 



We note that combined with the loading in a moving lat- 
tice described in section^ we can therefore prepare the 
atoms in any Bloch state \n,q). 

At the end of the acceleration period we ramp down 
the intensity of the lattice in 200 /is, while still moving 
at v^. After a 1.2 ms time-of- flight we take an absorption 
image of the cloud. A series of pictures corresponding to 
different final lattice velocities is shown in fig. 

Those pictures show that if the final lattice velocity 
remains within the first Brillouin zone (that is \vf\ < v^) 
the cloud comes back to rest in the laboratory frame after 
the adiabatic ramping down of the lattice. This behavior 
is now well understood in light of section IIIII On the 
other hand, each time the lattice final velocity reaches 
{2m-\-l)vr {m being an integer) the atom momentum, af- 
ter ramping-down the lattice, in the lab frame, increases 
by steps measured to be around 2Tik. This momentum re- 
mains constant for any lattice velocity between (2m+l)i'r 
and (2m + 3)i'r. 

As another way to understand this behavior in the first 
Brillouin zone, we again note that when the lattice, mov- 
ing with constant speed = -q^/M^ is ramped down 
adiabaticahy the velocity of the atoms with respect to 
the lattice varies from ^^(^f) to q^/M when the depth 
of the lattice goes to 0. The velocity of the atoms in the 
lab frame is thus g'f/M + 'Uf = 0. 

On the other hand, if the final velocity is, say, between 

and 3vr, the velocity of the atoms in the lattice frame 
is no longer q^/M after ramping down the lattice, but 
((7f + 2hk)/M. For example if the lattice is accelerated 
to Vi = 2.5vr, on ramping down, the velocity in the lat- 
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tice frame is —0.5vr. When the depth of the lattice ap- 
proaches 0, the velocity of the cloud in the lab frame thus 
goes to (gf + 2hk)/M — q^/M = -\-2vr. This explains the 
jump in momentum observed each time the velocity of 
the lattice reaches an odd number of recoil velocities. As 
in section^ the final momentum is independent of the 
intermediate lattice intensity. We have repeated the ex- 
periment for V = 1.5£^r, 5£^r and SEj- and found exactly 
the same behavior, apart from the small non-adiabaticity 
at the edge of a Brillouin zone. In the case of a shallow 
lattice, we interpret this jump in momentum in the labo- 
ratory frame as a first order Bragg diffraction: when the 
velocity of the lattice reaches an odd integer multiple of 
Vr thus statisfying the Bragg condition, the momentum in 
the lab frame changes by 2hk in the same direction as the 
acceleration. This Bragg diffraction is evidenced by the 
fact that the state of the atoms in the lattice connects 
to a different free parabola when the lattice is ramped 
down, as seen in fig. H 

One should not be misled by the fact that the conden- 
sate is back at rest in the laboratory frame when \v^\ < v^. 
The cloud has been displaced, dragged along by the lat- 
tice. The displacement is given by 



dE 
dq 



{q{t)) - q{t)/M dt 



(1) 



where r is the duration of the experiment (600 /is). 
For lattices deeper than about 3£^r, the derivative al- 
most vanishes and we approximate the displacement by 
X = (|raccei + Trampdown). In Order to determine 
whether the observed jump in momentum is exactly 2hk 
at the crossing of the edge of the Brillouin zone, one has 
to subtract this displacement due to the dragging of the 
lattice. This is shown in fig. The circles are the ac- 
tual positions of the center of the cloud in the lab frame. 
The crosses represent the positions corrected by the dis- 
placement due to the dragging. The dispersion of the 
data on a given plateau is due to a fiuctuation of the po- 
sition and velocity of the condensate, with RMS values 
of about, respectively, 10 jam and 0.03^'r. 

We next consider essentially the same experiment ex- 
cept that we now load the condensate in a lattice already 
moving with an initial velocity Vi = — Referring to 

figure^we see that adiabatic loading (100 /isec) prepares 
the atoms in the Bloch state \n = 1^ q = l.bhk). When 
the lattice is accelerated for 400 /isec in the positive di- 
rection in the lab frame, the atoms follow the first band 
and the quasimomentum in the lattice frame decreases 
linearly with time. Figure [5Ji shows the position of the 
cloud in the lab frame after the adiabatic ramp down of 
the lattice (100 /isec) and the subsequent 1.2 ms time of 
fiight. In figure [SJd we show the average momentum of 
the cloud. This includes compensation for the dragging 
of the atoms during the time the lattice is on (600 /isec), 
as described earlier. Figure [SId shows an alternation of 
—2hk and -\-4:hk momentum jumps in the lab frame. Ac- 
cording to the interpretation in terms of Bragg diffrac- 
tion, when the final velocity of the lattice reaches —v^ (or 
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FIG. 5: (Color) Acceleration of the lattice with atoms ini- 
tially loaded into the band n = 1, starting from q = 1.5hk. 
The 13Er lattice is raised adiabatically, accelerated, and then 
turned off adiabatically. Fig-I^i- shows the images of the cloud 
in the lab frame after a 1.2 ms time of flight for various final 
lattice velocities. Fig. \^ presents the momentum deduced 
from the positions of the cloud, after correction for the drag- 
ging. 



the quasimomentum reaches -\-hk)^ the atoms undergo a 
first order Bragg diffraction in the direction opposite to 
the acceleration of the lattice in the lab (equivalently 
they change from the Ohk to the -\-2fik free parabola, 
see figure 1). After being adiabatically released from the 
lattice, they now travel at —2vr in the lab frame, in the 
direction opposite to the acceleration. Despite the lattice 
being constantly accelerated in the direction of positive 
momentum, at this stage the atoms gain a momentum in 
a direction opposite to the acceleration. Further acceler- 
ation to Vf = leads to a second-order Bragg reflection 
that gives an impulse of +4fi/c, in the direction of the 
acceleration (corresponding to a change from the -\-2hk 
parabola to the —2hk parabola in the lattice frame, see 
fig.l). After adiabatic release the atoms' momentum in 
the lab frame is -\-2hk. 

As a conclusion of this section, we discuss the differ- 
ence between the experiments presented here and earlier 
experiments investigating Bloch oscillations. In refer- 
ence 0, for example, using the sudden turn-off method, 
the authors present the variation of the mean velocity 
of the atoms in the lab frame, after having accelerated 
the lattice. Their figure 2a shows steps of amplitude 2vr 
(note that their Vb = 2vr). The sharpness of the steps 
depends on the depth of the lattice and becomes more 
gentle when the lattice gets deeper (see their figure 2c) 
and figure 13 (upper panel) of g^]). On the other hand 
our figure^ exhibits sharp steps very similar to figure 
2a of p taken with a 0.29£^r deep lattice, despite the fact 
we were using a ISE^ deep lattice, the same as fig. 13 
of 3] . The adiabatic turn off method with any depth lat- 
tice thus produces results equivalent to the sudden turn 
off method using a vanishingly small lattice depth. This 
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FIG. 6: Same experiment as in figure2]and fig. IS) Velocity of 
the atoms in the lattice frame after acceleration of the lattice, 
for different final lattice velocities. The velocities are deduced 
from fig.EJ) and 13) by subtracting the velocity of the lattice. 
In fig. EJl the atoms are prepared in the ground state band, 
whereas they are prepared in the first band for the results 
of fig. ISJ). The plain lines are the theoretical group velocity 
calculated for a 0.1 Er deep lattice. 



is because when we turn the lattice intensity off adiabat- 
ically the states connect continuously to the Bloch states 
for a vanishingly shallow lattice. Figure IHt shows the ve- 
locity of the atoms with respect to the lattice when the 
lattice is off. This is equivalent to figure 2b of 5], with 
even sharper transitions. Our transitions are neverthe- 
less not infinitely sharp because we are not adiabatic very 
close to the zone boundary, as explained earlier. 

Based on this discussion and the data of figure we 
can infer what Bloch oscillations would look like in a 
weak lattice for a Bloch state in the first excited band. 
The velocity of the atoms in the lattice frame is presented 
in figure int). This figure was again obtained from figure 
ISJ) by subtracting the velocity of the lattice. Note that 
in contrast to the usual Bloch oscillations in the lowest 
band, here the Bloch oscillations in the first excited band 
consist of a series of first and second order Bragg diffrac- 
tions, at integer multiples of hk (half a reciprocal lattice 
vector) each reversing the velocity in the lattice frame. 
The first order Bragg diffraction changes the velocity in 
the direction of the force acting on the atoms in the lat- 
tice frame, whereas the second order Bragg diffraction 
changes the velocity by twice as much in the direction 
opposite to the force. This is in contrast with Bloch 
oscillations in the ground band where Bragg reflection 
occurs at multiples of 2hk (one reciprocal lattice vector), 
always in the direction opposite to the force. 



V. ACCELERATION OF ATOMS WITH A 
BROAD QUASIMOMENTUM DISTRIBUTION 

In a last set of experiments, we investigate the behavior 
of the atoms under acceleration of the lattice when the 
atoms do not occupy a single quasimomentum, but rather 
have a wide spread of quasimomenta. 

In order to prepare a broad distribution of quasimo- 
menta, we first reproduce the experiment of while 
the magnetic trap is still on at a relatively high mean 
oscillation frequency of 100 Hz in order to increase the 
interaction strength, we adiabatically turn on a 5Ej- deep 
lattice in 300 /is. We then suddenly turn off the magnetic 
trap 24] and let the atoms sit in the lattice for a duration 
ranging from 100 jas to 12 ms. We follow the evolution 
of the quasimomentum distribution of the atoms in the 
lattice by adiabatically turning off the lattice (in 300 /is) 
and taking an absorption image of the cloud after a 3 ms 
time of flight. The results are shown in fig. [7| [7^ shows 
an image of the undisturbed condensate after the time 
of flight, as well as the density profile along the lattice 
direction, integrated along the perpendicular direction, 
with no lattice having been applied; in[7|3, the lattice has 
been switched off suddenly, immediately after adiabatic 
loading of the lattice 25]. The momentum components 
at ±2hk appear and provide a calibration of the scale; 
in the two last pictures, [7t and [7|i, the condensate sits 
for respectively t = 0.5 ms and t = 9 ms in the lattice 
after which the lattice is ramped down in 300 /is. After 
9 ms of evolving time, when ramping down the lattice, 
the momentum distribution of the cloud looks essentially 
like a convolution of the profile of figure [7^ with an al- 
most uniformly populated momentum distribution with 
a 2hk width. This corresponds to an almost completely 
filled first Brillouin zone. More quantitatively, from the 
integrated profiles we calculate the rms width of the mo- 
mentum distribution of the atoms in the lab frame, after 
ramping down the lattice (see figure |H|). Since the ob- 
served distributions result from convolution of the quasi- 
momentum distribution with the distribution represented 
by figE^, we could in principle deconvolve them in order 
to get only the contribution of the quasimomentum. In- 
stead, as a reference, we show in fig. |Slthe expected rms 
width, convolving the experimental distribution of fig.[7tt 
with a quasimomentum distribution filling the first Bril- 
louin zone. 

An explanation for this broadening comes from the 
mean field inhomogeneity across the cloud. In the mag- 
netic trap, the chemical potential is independent of the 
position. When the lattice is super-imposed onto the 
magnetic trap, this is roughly still the case, provided that 
the lattice is not too deep. When the magnetic trap is 
suddenly turned off, the magnetic energy no longer com- 
pensates for the mean field energy and the chemical po- 
tential varies quadratically along the direction of the lat- 
tice. The rate of change of the phase difference between 
two neighboring sites then varies linearly along the lattice 
direction. This inhomogeneity of the density across the 
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FIG. 7: Dephasing of a condensate sitting in a 5Er deep lattice 
for a time t. The time of flight has the same 3 ms duration for 
all the images. The right column shows the density profiles 
of the images in the left column, integrated perpendicular to 
the axis of the lattice. In figlZK? no lattice was applied. In 
figEb the lattice was suddenly switched off. In[7|: and [31 the 
cloud stays in the lattice for 0.5 and 9 ms respectively. 
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FIG. 8: Evolution of the rms value of the observed the quasi- 
momentum distribution of the atoms as a function of the time 
spent by the atoms in the lattice. The dotted line repre- 
sents the expected width when we convolve the distribution 
of fig. |7Jl with a uniformly populated first Brillouin zone. 



condensate results in a different phase evolution at each 
lattice site and consequently in an effective dephasing 
of the single particle wave-function. Remembering that 
the quasimomentum characterizes the phase difference 
from one site to another, the apparent randomization of 
this phase leads to a broadening of the quasimomentum 



distribution. Roughly speaking, when the phase differ- 
ence between adjacent sites at the edge of the condensate 
reaches 27r, the wave function of an atom looks dephased, 
meaning that it is a superposition of all quasimomenta 
in the first Brillouin zone. This phase difference between 
neighboring sites at the edge of the condensate being on 
the order of jat/ (Nh) (where 2N-\-l is the total number of 
sites), after a time evolution of duration the time scale 
for dephasing is ~ Nh/ji. Numerically, the condensate 
had a chemical potential ja/h = 5 kHz. The estimated 
dephasing time (the time to create a 27r phase difference 
between adjacent wells) is then 8 ms, which is approx- 
imately equal to the observed time to fill the Brillouin 
zone. This treatment neglects tunnelling between lattice 
sites, which would tend to equalize the phases. However, 
we calculate a tunnelling rate of 27r x 1.5 kHz for a 5£^r 
deep lattice, that is, the well-to-well tunnelling rate is 
faster than the differential well-to-well phase evolution. 
Therefore, our simple picture of dephasing is question- 
able, although it seems to give a reasonable description 
of the experiment. We believe this point deserves further 
study. (A more detailed study of some aspects of mean- 
field dephasing in a lattice has been performed in |j3-) 

As an additional, albeit equivalent, demonstration for 
the randomization of the phase, we look for diffraction af- 
ter letting the condensate sit for a period of time. When 
we suddenly turn off the lattice, we do not observe re- 
solved diffraction peaks when the atoms have spent more 
than 2 ms in the lattice. We conclude that even though 
we may not have uniformly filled the Brillouin zone after 
2 ms, we broaden the quasimomentum distribution suf- 
ficiently that diffraction is not evident. We note that in 
reference 15j, the authors saw a diffraction pattern from 
an array of about 30 independent condensates. The dif- 
ference is in their smaller number of lattice sites, and may 
also be influenced by differences in experimental details 
such as the optical resolution for observing the diffrac- 
tion pattern, the number of diffraction peaks, and the 
fact that the diffraction of ref. appears not to be ob- 
served in the "far field" ^ . We apply the theory of pj| 
to our about 80 interfering condensates (assuming they 
are indeed independent which is only partially valid in 
our case) and found no diffraction pattern, in agreement 
with our observations. 

We finally turn to the behavior of the dephased cloud 
under acceleration. After letting the cloud sit in a bE^ 
deep lattice for 5 ms, more than a sufficient time to 
broaden the quasimomentum distribution enough that 
the diffraction is unresolved, we accelerate the lattice in 
500 /is to a chosen final velocity. After this acceleration 
period, we ramp down the lattice depth to zero in 100 fis 
and allow for a 3 ms time of flight, as described earlier 
in the paper. The results are presented in fig. |H1 In this 
figure, the mean momentum of the cloud after ramping 
down the lattice intensity shows no sign of the plateaus 
seen in fig. ^ This mean momentum is proportional to 
the lattice velocity, in contrast to the behavior described 
in section HVl In fact the atoms are dragged at the veloc- 
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FIG. 9: Acceleration of a cloud of incoherent atoms in a 5Er 
deep lattice. The circles represent the momentum of the cen- 
ter of mass of the cloud, in the lab frame, after acceleration 
of the lattice. The adiabatic ramp down is followed by a 3 ms 
time of flight. The line represents the velocity of the lattice. 



ity of the lattice, which means that in the frame of the 
lattice their motion is frozen. 

We reconcile this more intuitive behavior with the odd 
behavior of section IIVI by assuming that the first Bril- 
louin zone is completely filled, and considering a small 
component of the quasimomentum distribution, centered 
around go- Upon acceleration, this population does un- 
dergo Bragg diffraction when the velocity of the lattice 
reaches qo -\- hk^ and exhibits the same step behavior as 
the one seen in figure ^ the only difference being that 
the horizontal axis is shifted by qq. All the other quasi- 
momentum components are also Bragg reflected but at 
different velocities of the lattice. When one averages the 
different "staircase" patterns like the one shown in fig. ^ 
for all the quasimomenta in the first Brillouin zone, the 
average velocity in the lab frame is the velocity of the 
lattice. Another way to understand this is to calculate 
the average group velocity for a uniformly populated first 
Brillouin zone. That velocity is proportional to the aver- 
age of the slope of the Eo{q). As the band is symmetric 
with respect to g = 0, this average velocity with respect 
to the lattice is zero and there is no motion of the center 
of mass of the cloud with respect to the lattice. 

We now compare our above results with two recent 
experiments looking at thermal bosons and degenerate 
fermions in an optical lattice . In , a conden- 
sate surrounded by its thermal cloud is created in a sev- 
eral £^r-deep lattice and a magnetic trap. The center of 
the magnetic trap is then shifted and the subsequent be- 
havior of the two components is monitored. The authors 
observed that whereas the thermal component is pinned 
and does not oscillate in the magnetic trap, the conden- 
sate does oscillate with an oscillation frequency modified 
by the effective mass of the atoms in the lattice. The 



authors proposed an explanation based on the superflu- 
idity of the condensate that allows it to move through the 
corrugated potential created by the lattice, whereas the 
thermal cloud does not move due to its non-superfluid 
nature. In light of the experiment we described above 
in this section, we propose an alternate explanation for 
those results, based only on single-particle band structure 
theory. In the experiment of 16], the condensate is pre- 
pared directly in the lattice, and occupies the Bloch state 
|n = 0, g = 0). Let us now assume that the thermal com- 
ponent has a temperature that corresponds to an energy 
between the ground state band and the first band of the 
lattice. The ground band is then almost uniformly pop- 
ulated, meaning that the single-particle wave-function of 
the atoms effectively contains all quasimomenta in the 
first Brillouin zone. Shifting the position of the magnetic 
trap is equivalent to accelerating both the lattice and 
the trap with respect to the lab and therefore equivalent 
to applying a uniform force to the atoms. As described 
above, the thermal cloud filling the Brillouin zone does 
not move with respect to the lattice (figure 13 . This 
is what the authors of reference observed and it is 
completely consistent with a single-particle description, 
without any reference to superfluidity or critical velocity, 
phenomena dependent on interactions. 

We finally discuss briefly the recent experiment deal- 
ing with degenerate fermions in a one-dimensional optical 
lattice 17., In reference ^ Fermi sea of ^^K is 
produced in an optical lattice. The authors observe the 
absence of peaks in the diffraction pattern after sudden 
turn off of the lattice. This implies that the Fermi mo- 
mentum is comparable to or larger than hk so that the 
quasimomentum extends throughout the Brillouin zone, 
similar to our dephased cloud of Bosons. In the same 
work the authors repeat, with the Fermi gas, the experi- 
ment of reference 16j where they shift the magnetic trap 
with respect to the lattice. Consistent with our single- 
particle interpretation of the experiment (and with the 
single particle interpretation given in 17]), they do not 
observe oscillation of the Fermi cloud in the magnetic 
trap. In the authors again produce the Fermi sea 
in a lattice but this time they only partially fill the first 
Brillouin zone. As a result they do observe a diffraction 
pattern consisting of resolved peaks when suddenly re- 
leasing the atoms from the lattice. They also observe 
Bloch oscillations in their vertical lattice, due to gravity, 
as should be the case for a partially filled Brillouin zone. 
As ultra-cold indistinguishable fermions are essentially 
non- interacting (no 5-wave collisions), the experiments 
of 0, illustrate single particle, i.e. non-interacting 
particle, behavior of a cloud of cold atoms in a lattice, as 
those authors point out. Collisions imply a coupling be- 
tween quasimomentum states and thus the failure of the 
single particle (Bloch theory) description. The behavior 
of ultra-cold fermions is identical to our experiment and 
the experiment of |3 with interacting bosons, when the 
influence of interactions is negligible on the time scale of 
the experiment. It is particularly striking that Fermions 



10 



and Bosons can behave exactly in the same way under 
some circumstances: ah that matters is the way the Bril- 
louin zone is fihed, although the way this filling occurs 
may depend on the quantum statistics. 

VI. CONCLUSION 

In summary, we have presented a series of experiments 
in which a condensate is adiabatically loaded into an op- 
tical lattice, preparing the atoms in a single Bloch state. 
In a first set of experiments, the lattice is initially moving, 
and the atoms come back to rest in the lab frame after the 
adiabatic turn off of the lattice, leading to non-intuitive 
behavior for this "quantum conveyor belt". In a second 
set of experiments, we act on the prepared quasimomen- 
tum distribution by accelerating the lattice. We then 
analyze the new quasimomentum distribution by again 
adiabatically ramping down the lattice, and again ob- 
served non-intuitive behavior. We observe discrete jumps 
in the resulting momentum distribution, depending upon 
the velocity of the lattice. These jumps are reminiscent 
of Bragg diffraction at each avoiding crossing due to the 
laser coupling and are equivalent to Bloch oscillations. 
In a last set of experiments, we let the initial quasimo- 
mentum distribution evolve under the influence of inter- 
actions between the atoms, leading to the filling of the 
first Brillouin zone. The resulting cloud now exhibits a 
different behavior under acceleration of the lattice, i.e. 



the cloud appears to be frozen in the frame of the lattice. 
Finally we showed the similarities between the behavior 
of a cold thermal cloud and that of a cloud of degener- 
ate fermions in an accelerated optical lattice, when the 
quasimomentum extends throughout the first Brillouin 
zone. 

Among the issues that we believe deserve further stud- 
ies, both experimentally and theoretically, is the compe- 
tition between phase winding and tunnelling, that is to 
say how atoms lose their well- to- well phase coherence. 
Furthermore, we considered in this paper the dephas- 
ing of the wave-function due to the density profile of the 
cloud, but the quantum fluctuation of the atom number 
in each well is also a source of effective decoherence that 
should be explored. We also emphasize that the time 
scales of our experiments are very short with res2ect to 
other experiments, such as the one described in [13. Uq. 

Finally we note that the method of section HVl could 
be useful for precision measurements, as described in 
ref. 0. 
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